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A new class of nodal topological excitations in a two-dimensional Heisenberg model 
is studied. The solutions correspond to a nodal singular point of the gradient field of 
the azimuthal angle. An analytical solution found for the isotropic case. An effect of 
in-plane exchange anisotropy is studied numerically. It results in solutions which are 
analogues of the conventional out-of-plane solitons in the two-dimensional magnets. 
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The known topological solutions in two-dimensional (2D) Heisenberg model belong 
to classes of the homotopic groups isomorphic to the group of integers such as 7i 2 (S 2 ): 
Belavin-Polyakov (BP) solitons^ and easy-axis solitons 2 ; to the relative homotopic group 
7^2{S 2 , S 1 ): an out-of-plane (OP)^ and Takeno-Homma (TH)^ solitons; and to the group 
7Ti(5' 1 ): Kosterlitz-Thouless (KT) vortices. They correspond to a map of a spin order pa- 
rameter space onto a sphere S 2 homeomorphic to the 2D plane or a circle 5 1 — . 

In this paper we use an another way of a search of non-trivial topological excitations. By 
starting from the non-linear equations for (9, <^)-fields describing the dynamics of a classical 
2D isotropic Heisenberg magnet of a spin S with an exchange J 6 

= -KSsmO^ + JS 2 (a9 - cos 9 sinfl (w) 2 ) , (1) 
f)9 

= hS— + JS 2 (2 cos 9 (V0W) + sin 9A(f)) (2) 



dt 

one may see that in the stationary case the field variable 9 is determined by only the 
gradient field if) = V0. The Euclidean 2D-plane is homeomorphic to the sphere S 2 with one 
punctured point. According to the Hopf theorem^ Eulerian characteristic of triangulated 
surface equals to the sum of indices of singular points of the vector field on the surface. An 
Eulerian characteristic of the sphere is equal 2. In the simplest case one may suggest that 
one singular point with the Poincare index +1 is placed in the south pole and the other one 
at infinity with the same index corresponds to the punctured point in the north pole. The 
space configuration of the #-field will depend on the kind of singular point of if): a center, 
a node or a focus. The center singularity (Fig. la) corresponds to the solution of Eq.(J2J) 
A<p = and the particular solution = gtan -1 (j^j results in the well-known solitons with 
an axial symmetry listed above. One may expect that the focus singularity (Fig. lb) will 
correspond to a spiral spin arrangement. 

— # 

The investigation of the paper is devoted to the nodal point of the vector field if) (Fig. lc), 
where if) has a maximal (minimal) value. This means the choice of the following parametriza- 
tion (f) = cf) (r), i.e. the azimuthal angle changes along the radial direction in the plane. From 
the equation (J2J) one obtain 



dip q 



dr r sin 2 9 ' 

this determines the radial dependence of the <y9-field 

dr 



(3) 



/At 
, ■ 2 Q > ( 4 ) 
r sin 9 
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where the notation a is used for the lattice unit; ifo is an initial value. The equation for the 
9- angle 



A9 



cos 9 q 







sin 3 9 r 2 

may be integrated exactly and results in the scale invariant solution 



(5) 



9{r) 



cos 



\ 



\\Q\ 



±1. 



(6) 



This presents annuluses divergent logarithmically from the center (Fig. 2) and looks like 
a "target" with annular domains of magnetization. The parameter Q 2 > q 2 governs the 
amplitude of the oscillations, R is a scale factor, the sign p is a polarity of the solution. 
The continuum description is valid just for distances greater than the lattice unit a, the R 
value determines the boundary value 9q = 9\ r=a . As is seen from Eq. (JH} the 9 can not 
take the values and 7T. The plane chirality is determined by the sign of parameter q; the 
corresponding in-plane spin texture is presented in Fig. 3. 

This solution is a counterpart of the BP soliton because it has no definite localization 
radius and it is scale-invariant. In contrary to the BP soliton the energy of the found solution 

L r 



E 



JS 2 



V0 



q 



r 2 sin 9 



dr = nJS 2 Q 2 log I - 

.a 



(7) 



has no finite value and reveals the KT logarithmic behavior with an increase of the system 
size L. 

An availability of the pair of stationary solutions resembles in somewhat an exsistence of 
fundamental system of linear second-order differential equation. In the study of topological 
excitations in the classical XY model we deal the last situation: in Ref. 9 just one of the 
harmonic function = gtan -1 has been considered as physically reasonable; another 
solution = — q log(r) is exploited to obtain an effective interaction between vortices. 

As it follows from Eq.(|5|) the choice \Q\ — q realizes a pure in-plane arrangement 9 = | 
and insures a minimal energy E in the class of solutions. In the logarithmic scale x = log 
the parameter Q determines the wave length 5 — ^j, i.e. a distance between the nearest 
"crests" of the "target". The change of the angle tp on the scale S is A(f = qd that leads to 
the important relation for the small amplitude of the magnetization oscillations (q ~ Q) 



~2T 



g 

\Q\ 



(8) 
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i.e. the ratio qj \Q\ responds to the relative change of the azimuthal angle on the scale S. 

The criterion of topological stability of the BP soliton is rather simple^: there is an 
integer-valued topological invariant Q (degree of a map S 2 — * S 2 ) associated with the BP 
solution 9(f), 0(f) via 

Q = — I sm9(r)d9(f)d(f)(f). 

For the soliton with an axial symmetry = 0(a) (a is the angular polar variable) the 
topological invariant 

Q = ~ [cos 9(0) - cos#(oo)] = ±u, 

where v is a winding number. This describes a change of the angle at moving around the 
center r = 

(50 = 27TU 

and equals to degree of a map S 1 — ► S 1 . A non-zero density of an angular momenta 10 

L z = uhS(l — cos 9) 

and a zero value of the radial part of the momentum density P r = is a common point of 
the solitons with the center singularity: BP and easy-axis solitons, the OP and TH solitons, 
and also KT vortices. Thus, the winding number is associated uniquely with the angular 
momenta density of the magnetization field. 

For the solutions with a nodal singularity a situation is opposite: L z = and P r has a 
non-zero valued 

P r = TiS(l — cos 9)— = k -. 

dr r 1 + cos 9 

For a pure in-plane spin arrangement the q value determines a scaling factor 

log A = — , 
Q 

thus, that a change of the azimuthal angle on the scale A 

0(Ar) - 0(r) = 2tt. 

As is seen the q value is associated identically with the momentum density of the magneti- 
zation field. 

A special kind of nodal time-dependent solutions of two-dimensional Heisenberg model 
has been considered in Ref.— by means of the inverse scattering transform. Keeping time- 
dependent terms in the Eqs. (QQ |2J) one may obtain solutions with a finite energy and finite 
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localization radius. The radial behavior of the angle variables differs form the situation 
considered above, for example, the 9 reveals an exponential dependence at large distances 
in fixed time. The solitons are ring-shaped waves. Their localization radius and thickness 
grows with time linearly, whereas an amplitude is inversely proportional to the time. 

It is well known that an exchange anisotropy along z-axis J z > J± changes the asymptotic 
behavior at large distances from the power law decreasing to the exponential ones and results 
in the easy-axis solitons. An analogous effect occurs in the considered case but only for the 
plane exchange anisotropy J± > J z . The static equations for the case with an account of an 
external magnetic field along z-axis h may be written as 

= J ± (4cos0sin0 + cos 2 9A9) - J ± cos 9 sin 9((V9) 2 + (V0) 2 ) 

- 4 J z cos 9 sin 9 + J z cos 9 sin 9{V9f + J z sin 2 9A9 - - sin 9, (9) 
O = sin0A0 + 2cos0(V0V0). (10) 

The equation (jl(J|) is the same as in the isotropic case and results in the relation (@J). An 
analysis of asymptotic behavior at infinity yields the boundary value 

from which one obtain 9 = | at h — and there is no static solutions at h and 
J_i_ = J z . However, in the last case there is a dynamical solution ip = ip(r) + out with the 
Larmor resonance frequency. 

For a zero magnetic field h an asymptotic behavior at r — ► oo 

« = ! + » 

may be obtained from the equation 

The solution is the McDonald's function 59 = K v (jj, v = iq.\J~^f with the localization 
radius 

(i3) 

A series expansion of the solution © near the point r = a 

9{r) w 9(a) + 9' (a) r. (14) 
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A numerical calculation of the Eqs. (JHJ ITUj) is made by the shooting method and an example 
of the nodal OP soliton is presented in Fig. 4a. The phase portrait of possible behavior 
of the solutions presented in Fig. 5. An arbitrary choice of the derivative 0' (a) results in 
energetically unfavorable oscillating solutions above or below 9 = | (left and right limit 
cycles) and only an unique choice of the derivative gives the soliton. Like in the case of 
the conventional OP with an axial symmetry an energy has a logarithmic divergence with a 
growing of the soliton size. 

Finally, we point out that an inclusion of an external magnetic field h changes the asymp- 
totic behavior of the soliton at infinity from the exponential into the power decay 
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1 J ± cos 9 q 
66 ~ 4J ± -J 2 sin 5 fl r 2 ' { ' 

where 8q is given by This is a counterpart of conventional TH solitons (Fig. 4b). 

In conclusion, the class of the nodal stationary states in the 2D Heisenberg model is 
investigated. The scale-invariant solution, a counterpart of BP soliton, is found. An account 
of in-plane exchange anisotropy yields the analogues of out-of-plane and Takeno-Homma 
solitons, however, in contrary to the conventional solutions, those which are considered in 
the paper have no axial symmetry. 

This work was partly supported by the grant NREC-005 of US CRDF (Civilian Research 
& Development Foundation), INTAS grant (Project N 01-0654), by the grant "Russian 
Universities" ( UR.01. 01.005). 
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FIG. 1: The types of vector field singular points on the sphere: center (a), focus (b), node (c). 

FIG. 2: The space magnetization distribution in the "target". The logarithmic scale is used. The 
inset shows a radial section of the plot. 

FIG. 4: The 9(r) dependence in the nodal out-of-plane (a) and Takeno-Homma (b) soliton. 



FIG. 5: The phase trajectories of possible solutions of Eq. ©. 



FIG. 3: In-plane spin arragement of a nodal singularity. The dotted circle displays spin directions 
on the equal distances from the center. 



